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INTRODUCTION 


Purpose 

It  is  the  purpose  of  this  study  to  utilize  existing  theory  and  data  on 
supersonic  flows  similar  to  jet  engine  combustor  burn-through  flames,  along 
with  the  measured  flat  plate  impingement  pressure  and  temperatures  of 
Ref,  1  to  dete  rmine  inter  alia  the  distributions  of  velocity,  Mach  number, 
pressure  and  stagnation  temperature  (corrected  for  radiation  and  conduction 
errors)  along  the  flame  centerline. 

Background 

Considerable  background  on  the  problem  of  jet  engine  combustor  burn* 
through  failures  is  given  by  Rust?  who  gives  the  following  description  of  a 
burn- through: 

"A  burn- through  occurs  in  an  engine  when  the  hot  combustion 
gases  within  a  combustion  chamber  are  deflected  from  their 
normal  path  and  impinge  on  the  wall  of  the  combustion  chamber, 
thus  causing  a  hot  spot  on  the  wall.  The  heat  on  the  wall  weakens 
it,  and  the  high  pressure  in  the  chamber  causes  a  bursting  type 
failure  at  this  point.  The  result  is  a  severe  high  temperature, 
high  velocity  flame  escaping  from  the  hole  in  the  chamber.  " 

In  Ref,  1,  the  results  of  an  experimental  study  of  a  barn- through  flame'  on  a 
J47-GE-25  jet  engine  combustor  are  presented.  In  this  investigation  the 
flame  .vas  simulated  by  (1)  making  an  irregularly  shaped  hole  in  the  flame 
tube  which  fits  inside  the  burner  can,  (2)  inserting  a  baffle  plate  between  the 
flame  tube  and  the  burner  can- -creating  a  low  pressure  region  afr.  of  the 
baffle  plate- -causing  the  flame  to  impinge  on  the  burner- can  wall,  and 
(3)  making  a  1-in.  -diameter  hole  at  the  impingement  point,  allowing  the 
combustion  gases  to  exhaust  to  atmospheric  pressure. 

A  number  of  candidate  firewall*  materials  were  placed  in  the  flame  to 
determine  how  long  they  could  withstand  the  severe  heating  conditions.  In 
addition,  plate  impingement  pressure  and  stagnation  temperature  were  mea¬ 
sured  in  an  attempt  to  determine  the  pertinent  flame  properties 5  e.  g.  ,  temper' 
ature,  pressure,  velocity,  etc.  These  measurements  were  made  with  a  pitot 
pressure-measuring  hole  and  a  thermocouple  inserted  at  the  center  of  a  flat 
plate,  positioned  normal  to  the  flow.  Measurements  were  made  on  the  flame 
centerline  at  various  axial  d'stances  from  the  burn-through  hole.  Although 
the  above  measurements  gave  the  pressure  and  temperature  experienced  by  a 
firewall  material  in  the  flame,  additional  interpretation  of  the  data  is  required 
to  determine  local  flame  properties. 


The  firewall  is  used  to  protect  other  parts  of  the  aircraft  from  a 
burn- through  flame. 


In  the  discussion  part  of  this  report,  a  technique  for  interpretation  of  j 

the  pressure  data  in  terms  of  the  local  flame  centerline  properties  is  | 

described.  This  is  followed  by  a  description  of  a  method  for  computing  the  * 

convective  heat  transfer  coefficients  needed  as  input  data  for  the  thermocouple 
radiation  and  conduction  corrections.  Discussed  next  are  the  effects  of  mixing 
between  the  flame  and  ambient  air  on  the  centerline  properties.  Finally  recom¬ 
mendations  are  presented  regarding  instrumentation  and  data  interpretation 
for  future  test  programs  to  be  carried  out  at  NAFEC,  similar  in  nature  to  the 
program  described  in  Ref.  1.  (All  references  are  in  Appendix  C.) 


DISCtJSSIOM 


Interpretation  of  FI?*-  Plate  Pressure  Data 

If  the  ratio  of  stagnation  pressure  in  the  combustor  to  ambient  pressure 
is  greater  than  about  ’  ,  8  (for  7=1. 3),*  the  flow  will  be  sonic  at  the  burn- 
through  hole.  Table  1  (All  tables  are  in  Appendix  B.)  gives  the  properties  in 
the  combustor  and  at  the  hole,  and  indicates  that  =  1  for  all  power  settings, 
except  50  percent.  It  is  itad  that  foi  power  settings  2:  70  percent  the  pressure 
at  the  hole  is  greater  than  1  atmosphere  (atm);  thus,  the  burn-through  flame 
is  actually  a  high-temperature,  underexpanded  sonic  jet. 

Shock  Locations:  In  order  to  gain  confidence  in  the  subsequent  analysis  of 
flame  properties,  shock  locations  obtained  from  flame  photos  for  the  80  per¬ 
cent  and  90  percent  power  settings^  were  compared  with  data  available  in  the 
lltei  acure^’^  on  underexpanded  sonic  air  jets.  The  results  of  this  comparison 
are  given  in  Fig.  1.  (All  illustrations  are  in  Appendix  A.)  It  can  be  seen  that 
all  data  on  first  shock  locations  are  consistent,  and  a  smooth  curve  can  readily 
be  drav  n  through  all  the  points.  NAFEC  data^  on  the  oth'*r  shceV  locations  are 
seen  to  je  consistent  with  the  ARAP  data."  Figure  1  was  used  tu  determine 
shock  locations  for  power  settings  where  photos  were  not  available.^ 


^  It  is  assumed  that  y  =  1.3  for  the  burn-through  flame. 

t  A  photo  of  the  flame  for  the  70  percent  power  netting  was  not  clear  enough 
to  accurately  define  the  shock  locations. 

t  One  method  of  extrapolating  the  results  of  Fig.  1  to  higher  engine  pres¬ 
sure  ratios  is  to  (1)  use  the  results  of  Ref.  3  to  determine  the  location  of 
the  first  shock,  and  (2)  extrapolate  the  curves  for  shocks  2,  3  and  4  paral¬ 
lel  to  the  curve  for  shock  i. 
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Flat  Plate  Normal  Shock  Corrections;  In  a  supersonic  flow,  there  is  a 
difference  between  the  impingement  pressiire  on  a  flat  plate  and  the  pressure 
recorded  by  a  pitot  probe  at  the  same  location- -because  the  normal  shock 
will  be  farther  from  the  plate  than  from  the  probe  (see  Fig.  2).  Thus,  the 
probe  will  experience  a  pressure*  consistent  with  the  Mach  number  at  the 
probe  location  (Mpj.Q{,g)  while  the  plate  pressure  will  be  consistent  with  a  Mach 
number  upstream  of  the  plate  location  (Mpiate)*  This  effect  has  been  observed 
experimentally  by  Snedeker  and  Donaldson*  (see  Fig.  39  of  Ref.  2).  In  addition 
to  the  above  effect,  the  shock  in  front  of  the  plate  may  coalesce  wich  the  shock 
in  the  jet,  giving  rise  to  a  different  shock  structure  than  would  exist  if  no  plate 
were  present. 


Because  of  these  effects  the  first  step  in  interpreting  the  flat  plate 
pressure  data  was  to  "convert”  these  readings  to  equivalent  pitot  probe  pres¬ 
sures;  i.  e. .  the  pressure  downstream  of  a  normal  shock  at  the  local  Mach 
number.  This  was  done  by  utilizing  the  photos  in  Ref.  2,  which  show  shock 
detachment  distances  for  a  flat  plate  located  at  various  positions  in  a  jet.  By 
interpreting  the  pitot  probe  data  of  Ref.  2  in  terms  of  local  Mach  numbers,  the 
curve  shown  in  Fig.  3  was  constructed.  This  figure  was  used  in  the  following 
way  to  determine  equivalent  pitot  probe  pressures  in  the  burn-through  flame: 


1.  From  the  flat  plate  pressure ^  (Table  2)  and  free  stream  stagna¬ 
tion  pressure^  (i.e. ,  upstream  of  the  normal  shock)  a  free  stream  Mach  num.- 
ber  was  computea  (from  standard  normal  shock  tables'*). 


2.  From  Fig.  3,  the  shock  detachment  distance  was  determined. 
This  fixed  the  point  in  the  flow  having  the  free  stream  Mach  number  found  in  1 , 
and  theiefore  fixed  the  pitot  pressures^  at  that  point. 


From  the  above  procedure,  pitot  pressures  were  determined  as  a  function  of 
distance  from  the  burn-through  hole.  These  are  plotted  in  Fig.  4,  v/hich  will 
be  discussed  in  more  detail  in  the  following  section  (Mach  Number  and  Pres¬ 
sure  Distributions). 


It  is  noted  that  both  probe  and  plate  pressures  are  stagnation  pressures 


downstream  of  a  normal  shock,  but  at  Mach  numbers  ^^late' 

respectively. 


All  flat  plate  pressures  could  not  be  used  because  of  probable  interference 
with  the  jet  shocks  at  certain  locations.  Pressure  data  not  considered 
valid  are  noted  in  Table  2. 


t 
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These  computatTons  will  be  discussed  in  the  following  section  (Mach  Num¬ 
ber  and  Pressure  Distributions). 


It  is  noted  that  pitot  pressure  is  defined  as  the  stagnation  pressure  down¬ 
stream  of  a  normal  shock  in  front  of  a  pitot  probe. 
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Mach  Number  and  Pressure  Distributions;  The  following  procedure  was 
used  to  determine  the  Mach  number  and  pressure  distributions: 

1.  The  shock  locations*  were  determined  from  Fig.  1. 

2.  Between  the  burn-through  hole  and  the  first  shock,  the  pitot  pr  es¬ 
sures  and  Mach  numbers  were  determined  using  the  method  described  under 
Flat  Plate  Normal  Shock  Corrections,  with  the  free  stream  stagnation  pres¬ 
sure  equal  to  the  chamber  pressure. ^ 


3.  Conditions  immediately  downstream  of  the  first  shock  (and"  all 
succeeding  shocks)  were  determined  by  assuming  the  Mach  number  equals 
unity^  behind  the  shock(s). 


% 


't.  From  the  stagnation  pressure  downstream  of  the  shock  (from  3 
above)  and  the  pitot  pressure^  (corrected  flat  plate  pressures)  the  Mach  num¬ 
ber  distribution  between  the  first  and  second  shocks  was  determined  (similar 
to  step  2.) 

5.  Static  pressures  were  determined  from  Ref.  4,  knowing  the  loczd. 
Mach  numbers  and  stagnation  pressures,  assoming  V  =  1.3.§ 

6.  The  above  procedure  was  repeated  until  the  computed  stagnatlcrx 
pressure  loss  across  the  shocks  was  not  sufficient  to  account  for  !he  measured 
flat  plate  pressures.  The  stagnation  pressure  vs.  distance  curve  downstreamof 
this  point  was  approximated  as  described  in  the  following  sample  calculation. 


I 
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For  the  50  percent  and  60  percent  power  settings,  there  should  be  no 
shocks  in  the  burn-through  flame  (see  Table  1);  the  flow  is  completely 
subsonic. 


There  are  no  stagnation  pressure  losses  (along  the  centerline)before  the 
firsT  shock. 

This  is  a  reasonable  assumption  since  the  Mach  nuriber  behind  the  shock 
is  approximately  equal  to  or  less  than  the  jet  exit  Mach  number  (M^^). 
Since  the  observed  shocks  are  oblique,  the  downstream  Mach  numbers  can¬ 
not  be  less  than  unity.  Thus  M2^®)  =  1  was  chosen  to  interpret  these  data. 


Specific  heat  ratios  should  increase  towards  1.4  with  increasing  distances 
from  the  burn-through  hole  as  ambient  air  mixes  with  the  combustion 
gases.  However,  for  the  range  of  Mach  numbers  under  consideration,  the 
effect  upon  the  flame  properties  is  negligible. 
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Results  of  calculations  using  the  above  procedure  are  given  in 
Figs.  la-4h.  A  detailed  example  for  the  95  percent  power  setting  follows 


Before  First  Shock 

•  pQ  =  4.5  atm  (from  Table  1) 

.  at  X  =  1  in. ,  Popp  =  3.0  atm  (from  Table  2) 

P 

o 

FP 

.  From  Ref.  4,  for  V  =  1.3;  — p -  =  0.  67,  gives  M  =  2.05 


g 


•  From  Fig.  3,  A/  d  =  0.  25 

•  Mjjj  =  2.05  is  then  located  at  x  =  0.75  in.  (see  Fig.  4a) 

•  From  Ref.  4,  p/Pq  =  0.  12,  thus  p  =  0.54  atm 

•  The  Mach  number  curve  is  faired  through  three  points; 

^ex  =  1  at  X  =  0,  and  2  points  (x  =  0.?5  and  1.2  in.)  deter¬ 
mined  in  the  above  r.iacjier. 


At  First  Shock 


•  From  -  2,3  (from  the  faired  Mach  Number  curve)  and 

=1.0,  the  stagnation  pressure  downstream  of  the 
oblique  shock  is  determined  (from  Ref.  5);  =  3.8  atiii. 

•  From  Ref.  5,  a  shock  angle  of  58*  is  found  for  this  pressure 

loss.  The  follow’ug  table  gives  a  comparison  of  first -shock 
ang;'  s  determined  by  the  above  method  anf’,  those  measured 
from  the  nhotos  o^  Ref.  1.  s 


K  fc 


Engine  Power  Setting,  Percent  RPM 

90 

30 


Shock  Angle,  deg. 
Predicted  Measured 


60 

64 


52 

58 


This  good  agreement  between  the  predicted  and  measured  shock  angles 
gives  additional  confidence  in  the  assumption  of  M2^°)  =  1.0. 
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Before  Second  Shock 


•  pQ  =  2.8  atm  (from  above) 

•  At  X  =  3  in. ,  ^opp  =  2.5  atm  (from  Table  2) 

^opp 

•  From  Ref.  4»  for  Y  =  1.3;  ~ —  =  0.66,  gives  M®  =  2. 07 

•  From  Fig.  3,  A/d  =  0.  25 

•  M®  =  2.07  is  then  located  at  x  =  2.75  in. 


This  procedure  is  continued  unt;J  the  calculated  stagnation  pressure 
loss  is  not  sufficient  to  account  for  the  measured  flat  plate  pressure.  *  At  the 
95  percent  power  setting  this  condition  occurs  between  the  third  and  fourth 
shock.  The  stagnation  pressure  curve  (Fig.  4a)  is  then  continued  by  connect¬ 
ing  the  pitot  pressures  behind  the  4th  and  5th  shocks  (since  M2'®'  =  1,  the 
stagnation  pressure  must  equal  the  pitot  pressure)  and  extrapolating  the  curve 
downstream  of  the  5th  shock. 


Velocity  Distributions:  Before  local  velocities  can  be  computed  the 
thermocouple  readings  must  be  corrected  for  radiation  and  conduction  losses: 
the«e  corrections  were  made  (see  Thermocouple  Corrections  for  details),  and 
the  true  gas  stagnation  temperatures  are  given  in  Fig.  4.  The  velocity  is  then 
computed  from: 


u  =  Ma 


=  M[7R(i)T„] 


where  T/Tq  is  a  function  of  Mach  number.* 


This  is  probably  due  to  an  increased  stagnation  pressure  loss  due  to  mix¬ 
ing  bebveen  the  ambient  air  and  burn -through  flame  (see  Influence  of  Mix- 
jjag  on  Flame  Properties).  Because  of  the  approximations  in  the  present 
analysis  and  the  inherent  inaccuracies  in  the  measured  flat  plat  impinge¬ 
ment  pressures,  this  point  cannot  be  determined  with  great  precision. 


Thermocouple  Corrections 

Radiation:  Because  of  radiation  from  the  thermocouple  junction  (and 
wires)  to  the  surroundings,  the  equilibrium  temperature  of  the  junction, 
is  less  than  the  true  gas  stagnation  temperature,  Tq.  However,  Tq  cam  readily 
be  determined  by  performing  a  heat  balance  (between  convection  to  the  junc¬ 
tion  and  radiation  from  the  junction).  Assuming  conduction  along  the  wires  is 
negligible*  we  get, 

h(To  -  T^)  =  0-£  (2) 

Eq.  (2)  is  used  to  solve  for  Tq,  once  the  heat  transfer  coefficiear,  h, 
is  known. 


Heat  Transfer  Coefficient  -  There  are  no  experimental  data  in  the 
literature  on  heat  transfer  coefficients  for  hot,  turbulent  supersonic  jets 
impinging  on  a  flat  plate.  However,  data  are  available  for  cold  subsonic  tur¬ 
bulent  jets  impinging  on  moderately  heated  flat  plates,^’’  and  a  (well -confirmed) 
theoretical  expression  is  available  for  stagnation  point  heat  transfer  in  laminar 
flow.®  These  pieces  jf  information  were  combined  in  this  study  to  obtain  the 
local  heat  transfer  coefficients. 

Doraldson,  et  al.,®  compared  their  measured  stagnation  point 
heat  transfer  coefficients  (for  turbulent  jets)  with  those  computed  from  laminar 
stagnation  poin,  theory®  and  showed  that  the  ratio  of  the  measured  to  calculated 
values  is  a  furxtion  of  distance  from  the  jet  exit  (and  not  a  function  of  Reynolds 
number).  By  extrapolating  the  data  of  Ref.  6. it  was  fj’ond  that  this  ratio  varied 
from  about  1,1  to  1,5  for  the  range  of  x/d  of  interest  in  the  present  study.  Thus 

the  laminar  stagnation  point  heat  transfer  coefficients  were  computed  from  the 
theory  of  Ref.  8,  and  all  values  were  multiplied  by  an  (average)  factor  of  1»3 
to  get  the  corresponding  values  of  h  for  turbulent  flowi. 

Values  of  h  for  laminar  flow  were  computed  from,® 

1/ 

t_  1  1/  1/  Pj  /J._v  "1 


Nu  =  ^  =  0.66  Re 

k 


[_u  Vdy /o  _ 


The  normalized  velocity  gradient;  i.  e. ,  the  term  in  square  brackets,  was 
measured  by  Snedeker  and  Donaldson^  (for  latios  of  jet  exit  to  ambient  pres¬ 
sure  from  1.0  to  3.5)  to  be  about  unity.  It  is  important  to  note  that,  for 


This  is  shown  under  Conduction. 

Use  of  an  average  factor  is  consistent  with  the  accuracy  of  this  technique 
for  calculating  turbulent  heat  transfer  coefficients. 
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readings  in  Table  3  the  actual  gas  stagnation  temperatures  were  computed; 
these  are  listed  in  Table  4  and  plotted  in  Fig.  4. 

Smooth  curves  were  drawn  through  the  data  in  Fig.  4  and  the 
faired  values  near  the  burn-through  hole  were  used  to  estin^ate  stagnation 
temperatures  in  the  combustor  for  all  power  settings;  these  are  listed  in 
Table  1.  Rust^  reports  that  information  furnished  by  the  J-47  engine  manu¬ 
facturer  (General  Electric)  indicates  a  maximum  stagnation  temperature  in 
the  combustor  of  3500*R.  In  Table  1  it  can  be  noted  that  this  is  exactly  the 
stagnation  temperature  estimated  for  the  95  percent  power  setting.  This  good 
agreement  implies  that  the  thermocouple  corrections  are  reasonable. 

Conduction:  Conduction  from  the  thermocouple  junction  along  the  wires 
can  influence  the  heat  balance  equation,  Eq.  (2),  if  the  conductive  heat  trans¬ 
fer;  i.  e. , 

‘Icond  “  ^  )j 

is  not  negligible  compared  to  the  convective  heat  transfer  to  the  junction, 

^cony  “  ^jC^o  “ 

Pertinent  data  on  the  thermocouple  are: 

•  Material:  Pt/Pt-13%  Rh  (Ref.  1) 

•  length  to  connectois  =  10  in.  (measured) 

•  d.„  =  0.  02in.  (measured) 

♦V  ' 

•  junction  (bead)  diameter  =  0.062  in.  (measured) 

As  an  example,  the  conductive  heat  transfer  for  a  thermocouple  reading  of 
3000®R,  and  an  (assumed)  temperature  gradient*  (dt/ds)-  =  300*R/in. ;  i,  e,, 

J 


This  assumed  gradient  is  conservatively  large  since  the  junction 
and  exposed  wires  are  both  subjected  to  the  convective  heating.  This 
would  make  the  temperature  gradient  at  the  junction  much  smaller 
than,  in  this  case,  300*R/in. ,  thereby  reducing  the  conductive  heat  lost 
even  further. 


a  linear  gradient  ever  the  10  in.  length  is 

^cond  ""  0*39  Btu/hr  (Eq.  4) 

For  Tq  =  3500*R,  =  3000*R  and  h  =  180  Btu/hr-ft’-*R,  the  convective  heat 

transfer  is 

•Iconv  =  Btu/hr  (Eq.  5) 

Thus,  even  using  a  conservatively  large  value  of  (dT/ds)j,  q^ond  case) 

is  only  20  percent  of  q^Q^v*  conduction  heat  loss  was  found  to  be  negligible 

for  all  test  conditions,  primarily  due  to  a  small  wire  diameter  and  relatively 
small  temperature  gradients  along  the  wire. 

Mass  Flow  Rates 

Following  determination  of  stagnation  temperatures  in  the  combustor,  the 
mass  flow  rates  can  be  coinputed  via, 

~  (Pex/^^ex)  “ex^ex 

Using  the  data  given  in  Table  1>  along  with  Eq.  (6),  gives  the  mass  flow  rates 
shown  in  Table  5, 


Influence  of  Mixing  on  Flame  Properties 

Mixing  between  the  burn -through  flame  and  ambient  air  is  the  most 
important  mechanism  by  which  stagnation  enthalpies  in  the  flame  decrease. 
Thus,  a  good  measure  of  the  effect  of  nuxing  on  the  flame  properties  along 
the  centerline  is  the  rate  of  decrease  of  stagnation  enthalpy.  *  Nondimensional 
stagnation  enthalpies  (the  usual  manner  of  presenting  high-temperature  mixing 
data)  are  plotted  in  Fig.  5  for  all  power  settings.  An  important  parameter  in 
the  present  study  is  the  distance  from  the  burn-through  hole  at  which  the  effects 
of  mixing  are  experienced  on  the  flame  centerline  (the  designated  core  region). 
The  length  of  the  core  region,  L,  is  arbitrarily  defined  as  the  value  of  x/d  at 


Enthalpy  is  csed  us  a.  mixing  parameter  rather  them  temperature  because 
the  specific  heat  of  the  flame  varies  over  the  temperature  range  of 
interest. 
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which.  -Hamb)  ~  0*95.  These  core  lengths  are  given  in  Table  6, 

together  with  the  angle  of  spread  of  the  mixing  region,*  i.  e.  t3-.n“^  d/ZL. 

In  order  to  extrapolate  these  test  results  to  other  systems  the  core  lengths 
are  correlated  with  enthalpy  at  the  burn -through  hole,  i,  e. , 

^  ~  Pex^^ex(^Oj.  "  ^amb) 

.  for  all  power  settings.  (See  Fig.  6).  The  data  can  be  fitted  reasonably  well 
by  a  straight  line  if  we  neglect  the  point  for  the  50  percent  power  setting. 

There  are  notmuch  data  available  on  high -temperature  turbulent  jets  to 
compare  with  the  results  of  the  NAFEC  tests.  However,  O'Connor,  etal. 
have  measured  the  axial  variation  of  stagnation  enthalpyt  in  turbulent  jets  of 
partially  dissociated  N2,  with  an  enthalpy  flux  of  about  5  X  10*  Btu/ft^-sec. 
These  data  are  shown  in  Fig.  5.  We  observe  that  the  corresponding  point  on 
the  core  length  vs.  enthalpy  flux  curve  (Fig.  6)  is  somewhat  lower  than  the 
NAFEC  data.  This  agreement  is,  howeven  good  considering  the  uncertainties 
in  estimating  core  lengths  from  the  enthalpy  data.  The  straight  line  correla¬ 
tion  should  be  useful  for  determining  core  lengths  and  mixing  region  angles 
(cf.  Table  6)  for  larger  burn-through  boles  with  greater  mass  flows. 

However  it  is  unadvisable  to  use  Fig.  6  to  determine  core  lengths  for  high 
compression  ratio  (say,  12:1)  engines,  for  which  the  burn-through  jet  bound¬ 
aries  are  highly  expanded.  Because  of  this  expansion,  the  core  lengths  should 
be  greater  than  shown  in  Fig.  6,  for  a  given  enthalpy  flux.  Interpretation  of 
highly  underexpanded  rocket  plume  data  would  be  needed  to  get  an  equivalent 
correlation. 

Some  additional  data^  on  stagnation  temperature  and  enthalpy  axial 
distributions  in  high  temperature  jets  are  given  by  Harsha.*^ 


*  Camion  must  be  exercised  in  attempting  to  compare  the  results  of  Fig.  5 
and  Table  6  with  corresponding  data  on  high  temperature  jets  issuing 
from  a  perfectly  expanded,  smoothly  converging  nozzle.  It  has  been  shown 
by  Donaldson  and  Gray^®  that  the  mixing  rates  are  higher  for  jets  emanating 
from  a  flat  surface  (which  would  be  similar  to  combustor  burn-throughs) 
than  for  the  above-mentioned  jets.  This  is  due  to  the  difference  in  air 
entrainment  patterns. 

t  These  data  are  for  a  jet  (at  a  pressure  of  1  atm)  emanating  from  a  flat 
plate,  so  they  should  compare  reasonably  well  with  the  NAFEC  data. 

J  There  was  insufficient  time  in  the  present  study  to  interpret  these  data. 
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Suggested  Experimental  Methods  for  Future  Test  Programs 

It  is  suggested  that  future  test  programs  on  higher  compression  ratio 
engines  be  designer,  to  obtain  data  on  local  burn-through  flame  properties 
(temperature,  pressure,  etc. )  and  local  heat  transfer  coefficients  (for  a 
geometry  similar  to  that  of  a  firewall  material  in  the  flame).  In  addition, 
analytical  models  which  can  be  used  to  predict  burn-through  flame  properties 
and  heat  transfer  coefficients  should  be  developed.  After  it  has  been  denAon- 
strated  that  predictions  and  test  data  are  in  good  agreement,  the  models  should 
be  used  to  develop  a  set  of  design  charts — giving  heat  transfer  coefficients, 
stagnation  temperatures  and  presstires  as  a  function  of  axial  and  radial 
distance  from  the  burn-through  hole,  for  parametric  variations  of  combustor 
stagnation  pressure,  temperature  and  hole  sizes.  These  are  the  necessary 
input  data  for  determining  the  response  of  firewall  imterials  to  these 
burn- through  flames.  B;  is  important  to  note  that  the  development  of  accurate 
analytical  models  will  avoid  the  necessity  of  obtaining  test  data  each  time  a 
new  combustor  is  designed  or  an  existing  system  is  re-designed.  These 
models  are  also  necessary  in  order  to  design  a  standard  flame  (which 
simulates  a  combustor  burn-through  flame)  to  test  firewall  materials. 

Specific  measurements  that  should  be  taken  in  future  test  programs 
are  listed  below. 

Burn- Through  Flame  Properties;  (These  shovild  be  taken  as  a  function 
of  axial  and  radial  distance  from  the  burn-through  hole  at  a  sufficient  nximber 
of  points  to  define  the  properties  within  each  shock  cell) 

1,  Pressure  -  Both  pitot  and  static  pressure  tubes  shotild  be 
used  to  record  both  the  stagnation  pressure  behind  a  normal  shock  and  the 
local  static  pressure  (cf.  Refs,  13  and  14).  These  data  are  then  used  to 
determine  local  Mach  numbers. 

2,  Stagnation  Temperature  -  A  stagnation  point  thermocouple 
should  be  used,  with  provisions  made  to  minimize  heat  conduction  errors 
(cf.  Ref,  15).  This  can  be  accomplished  by  using  very  small  lead  wires  and 
inserting  the  jtmction  far  enough  into  the  flame  so  that  the  junction  and  wires 
are  at  about  the  same  temperature, 

Hiotographs:  These  are  especially  important  for  determining  the  flame 
boundaries  and  shock  structure- -for  comparison  with  analytical  predictions. 
Color  photos  are  useful  for  estimating  the  extent  of  the  mixing  region— but 
black  and  white  photos  are  sufficient  for  defining  the  shock  structure. 
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Sampling;  It  would  be  veiy  useful  to  know  the  composition  of  the  burn- 
through  flame — particularly  to  determine  if  any  unburned  fuel  escapes  the 
combustor.  The  fiiel  could  (subsequently)  react  with  ambient  air  giving 
higher  flame  temperatures  than  would  otherwise  be  expected.  Gas  samples 
should  be  taken  via  a  sampling  probe  throughout  the  burn-through  jet  (in  the 
core  and  mixing  regions).  Samples  taken  at  the  burn- through  hole  will 
indicate  how  much  unburned  fuel  is  available  to  react  with  ambient  air.  A 
detailed  discussion  of  the  problems  involved  in  gea  sampling  and  the 
considerations  involved  in  designing  sampling  probes  are  given  in  Ret  16. 

Effects  of  Noncircular  Holes;  Since  a  burn- through  hole  will  jirobably 
not  be  circular  in  practice,  it  is  important  to  test  the  influence  of,  fer 
example,  an  elliptical  hole  with  jagged  edges,  on  flame  properties.  Signifi¬ 
cant  differences  between  test  results  on  circular  and  noncircular  holes  must 
then  be  accounted  for  in  predicting  actual  burn-through  flam'j  proT:;erties. 

Heat  Transfer  Coefficients;  These  can  be  determined  by  measuring 
local  heat  transfer  rates  and  wall  temperatures  fox  a  flat  plate  (simulating  a 
firewall  material)  normal  to  the  burn-through  flu  je,  (The  plate  impinge¬ 
ment  pressures  should  also  be  recorded  as  a  function  of  radial  distance  from 
the  plate  stagnation  point  since  these  are  the  actual  pressures  experienced  by 
firewalls. )  Two  possible  nonsteady-state  meth.'>ds  of  measuring  heat  transfer 
coefficients  are; 

1.  Calorimeter  Discs  (see,  e.  g.  Re£  6)  -  In  this  technique  the 
temperature-time  response  of  a  uniform  temperature  disc 
is  recorded.  The  heat  transfer  is  computed  from. 


q 


conv 


(8) 


then 

q 

conv 

^  T^  -  T  (9) 

o  w  '  ' 

where  T^  is  the  measured  3 'Agnation  temperature. 


2.  One-Dimensional  Plug  -  A  cylindrical  plug  is  inserted  in  the 
flat  plte,  with  one  end  exposed  to  the  gas.  The  surface  is 
insulated  over  its  length  and  temperature  of  the  exposed  end 
is  measured  as  a  function  of  time.  From  the  solution  to  the 
heat  conduction  equation  for  a  one-dimensional  semi-infinite 
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rod*  with  a  convective  boundary  condition,  and  'the  measured 
temperature  of  the  exposed  end,  the  heat  transfer  coefficient 
cam  be  determined  from. 


T  -  T.  r  ,  , 

To  -  Ti  a 


“*]  it^) 


Additional  information  on  this  technique  is  given  in  Ref.  18. 


*  The  plug  acts  as  a  semi-infinite  rod  because  data  are  taken  before  the 
temperature  of  the  back  end  can  increase. 


CONCLITSION 

Based  on  the  results  of  this  stud-/  program,  it  is  concluded  that: 

The  computed  pressures,  t'':mperatures  and  mass  flows  reported 
herein  closely  define  the  environment  experienced  by  any  material  encounter¬ 
ing  a  combustor  burn- through  flame  froi:!  an  engine  operating  within  the 
pressure  range  specified  in  this  report. 
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OF  BURN-THROUGH  FLAME 


50%  power  setting 


VELOCITY,  ft/sec 
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NAFEC  DATA 


Correlation  developed  for  ratios 
of  jet  exit  to  ambient  pressures  £  2.5 


TABLE  i  -PROPERTIES  IN  COMBUSTOR 


Engine  Power 
Setting 
(%  RPM) 

Stagnation 

Pressure^ 

(atm) 

Stagnation 

Temperature^ 

(•R) 

Static  Pressure 
at  Hole^ 

(atm) 

Velocity 
at  Hole® 
(ft/sec) 

50 

1.5 

2300 

1.0*^ 

1800*^ 

60 

1.8 

2600 

1.0 

2400 

70 

2.4 

3000 

1.3 

2600 

75 

2.8 

3100 

1.5 

2600 

80 

3.2 

3200 

1.8 

2700 

85 

3.6 

3300 

2.0 

2700 

90 

4.1 

3400 

2.2 

2800 

95 

4.5 

3500 

2.5 

2800 

a  Taken  from  Ref.  1. 

b  Interpreted  from  thermocouple  data  corrected  for  radiation  and  conduc 
tion  losses  (see  Fig.  4). 

c  From  isentropic  flow  tables  (Ref.  4);  Y  =1.3. 


d 

e 


Mgx  ~  0.8.  =1.0  for  all  other  power  settings. 

From  Eq.  (la)  using  Y  -1.3  and  R  =  53.3  ft-lb£/lbjjj-®R. 


Preceding  page  blank 
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TABLE  5-  MASS  FLOW  RATES 


Engine  Power  Setting 
(%  RPM) 

Mass  Flow  Rate 
(Ib^j/sec) 

50 

0. 18 

60 

0.23 

70 

0.28 

75 

0.31 

80 

0.38 

85 

0.42 

90 

0.45 

95 

0.49 

a 


FromEq.  (6),  using  Y=  1.3  and 
R  =  53.  3  a-lbf/lb^-^R 
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